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Abstract 

In this paper, we present a mean field game to model the production behaviors 
of a very large number of producers, whose carbon emissions are regulated by gov¬ 
ernment. Especially, an emission permits trading scheme is considered in our model, 
in which each enterprise can trade its own permits flexibly. By means of the mean 
field equilibrium, we obtain a Hamilton-Jacobi-Bellman (HJB) equation coupled with 
a Kolmogorov equation, which are satisfied by the adjoint state and the density of pro¬ 
ducers (agents), respectively. Then, we propose a so-called fitted finite volume method 
to solve the HJB equation and the Kolmogorov equation. The efficiency and the use¬ 
fulness of this method are illustrated by the numerical experiments. Under different 
conditions, the equilibrium states as well as the effects of the emission permits price 
are examined, which demonstrates that the emission permits trading scheme influences 
the producers’ behaviors, that is, more populations would like to choose a lower rather 
than a higher emission level when the emission permits are expensive. 

Keywords. Mean field games. Producers’ behaviors. Emission permits trading. Fitted finite 
volume method. Equilibrium states. 

*This project was supported in part by the National Basic Research Program (2012CB955804), the Major 
Research Plan of the National Natural Science Foundation of China (91430108), the National Natural Science 
Foundation of China (11171251), the Russian Foundation for Basic Research (14-07-00075), and the Major 
Program of Tianjin University of Finance and Economics (ZD1302). 


1 



1 Introduction 


Due to the rapid development of economy and society, every one should analyse and get 
a clear understanding of a very complex system before making a decision. In particle physics, 
mean held theory can be regarded as a useful and effective tool to study the performance 
of large and complex stochastic model. It focuses mainly on a single particle and assumes 
that the interactions of this particle with its neighboring particles is determined by the mean 
held, in which the inter-particle interactions must be sufficiently “weak” or “regular” and 
each particle tends to be inhnitesimal, i.e., the number of particles tends to inhnity. 

Based on the above features of mean held theory, Lasry and Lions hrst dehned and 
developed mean held game in their three papers [1-3]. Diherent from the standard game 
theory in which the number of players is hnite, such as [4-7], mean held game studies the 
limit case of a game with N players as N goes to inhnity, which implies a continuum of 
agents. Owing to this, mean held game can be used to deal with the problems, which should 
be summed up by an untractable system of HJB equations in general diherential games with 
N players, where N is very large. 

Several researchers have studied the theory and applications of mean held game in recent 
years [8-12]. Especially, some interesting and meaningful works relating the mean held game 
to economics have been done. Gomes et ah present ehective numerical methods for two-state 
mean held games and discuss a number of illustrative examples in socio-economic sciences in 
[13]. Besides, in [14] Lachapelle et ah consider a technology choice problem with externalities 
and economy of scale by using a mean held game stylized model. They introduce a monotonic 
algorithm to hnd the mean held equilibria and describe the multiplicity of equilibria. 

However, to our best knowledge, there are very few studies on mean held games to take 
emission permits trading into consideration. For the past few years, the issues on climate 
change and emission reduction have attracted the attention of politicians and scholars from 
all over the world. In order to mitigate climate change and improve global environment, 
some emission permits trading markets have emerged and are developing prosperously. At 
the same time, a large number of articles have studied the emission permits price theoretically 
and empirically [15-19]. Moreover, some diherential game models about production decision 
also include the emission permits trading scheme [20, 21]. 

Motivated by the those mentioned above, in this paper, we build a mean held game 
model, in which the revenues from emission permits trading are included, to study the 
productive behaviors of a continuum of agents under the background of climate change and 
emission reduction. In our model, each agent is anonymous and the interactions among 
agents are mean held type. They can obtain an initial quota from the emission permits 
trading scheme, and purchase the emission permits from the market compulsively if the 
quota is insufficient or sell the unused emission permits to others. The equilibrium of mean 
held game can be reached by solving two coupled partial diherential equations, one of which 
is a backward HJB equation satished by the adjoint state, and the other one is a forward 
Kolmogorov equation satished by the density of agents. 

Some discussions about the numerical algorithms of the above coupled system have 
been made for the past few years [13, 14, 22-24]. All of these methods are based on a hnite 
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difference scheme. In this paper, we propose a so-called htted hnite volume method to solve 
the coupled equations model established by ourselves. The innovation of this method is in 
that it couples a hnite volume formulation with a htted local approximation to the solution. 
On the one hand, we implement the local approximation through solving a sequence of two- 
point boundary value problems dehned on each element. On the other hand, the hnite volume 
method possesses a special feature of the local conservativity of the numerical hux. The main 
advantage of this discrete method is that the system matrix of the resulted discrete equation 
is an M-matrix, which guarantees that the discretization is monotonic and the discrete 
maximum principle is satished. The hnite volume method, except for the fundamental huid 
dynamics problem in which it performs very well, has been used in many helds and is 
becoming more and more popular. See, for instance, [25] for degenerate parabolic problems, 
[26] for hyperbolic problems, and [27] for elliptic problems. 

The paper is organized as follows. In Section 2, a mean held game model is established, 
and the coupled partial diherential equations are presented. Then, a so-called htted hnite 
volume method is proposed for the discretization of the equations in Section 3. In Section 4, 
some numerical experiments are performed to illustrate the efficiency and usefulness of the 
numerical method, and the ehects of the parameters on the density of population are also 
showed in this section. Finally, concluding remarks are given in Section 5. 


2 The mean field game model 

2.1 The states of agents 

For the purpose of illustrating the interactions among the players in a commitment 
period, we propose a hnite-horizon mean held game framework. We focus on a very large 
economy and consider a continuum of agents. Each agent is a producer whose carbon 
emission is regulated by governments. In our settings, all the producers have the same 
capacity in production. They are anonymous but diherent in their initial production states 
which follow a given initial probability density function. In addition, the interactions among 
the agents are mean held type. That is to say, a given agent cannot inhuence the distribution 
of all players’ states and therefore the decisions of others by itself. However, it can produce 
an ehect on the information which is used by others to make decisions. These ideas, as well 
as the name of mean field, come from particle physics, and the particles are replaced by 
rational agents here. 

Let Q{t) denote the production of each agent during the period of [0, T], where T is the 
maturity of the game. This production leads to an amount of by-products, namely emissions 
E{t). Generally speaking, an increase in production can result in more emissions and vice 
versa. So, it is reasonable to use the emission as a state variable instead of production in this 
paper. In fact, the emission level can be also treated as a product portfolio. The dynamics 
of the agent’s emission is given by the following controlled process: 

dEt =-T{t,E)dt + (TdW+ dNt{Et), (1) 

where r(t, E) is a control variable, and can be interpreted as the level of emission reduction. 
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The term adW stands for the stochastic disturbance of emission resulting from the techno¬ 
logical innovations, market fluctuations, and some other uncertain factors. The constant a 
is a noise parameter and denotes the volatility of emission, and dW is the increment to the 
standard Brownian process. In addition, the emission Et is restricted in -^max] by fhe 

reflection part dNt{Et), in which Nt is a monotonic continuous nondecrease function. For 
more details about this formulation, please see [14] and [28]. 

Taking advantage of the dynamics of emission (1), we can obtain the forward Kol¬ 
mogorov equation satisfied by the density of agents m(t, E) for t G [0, T] and E G [i?min, T^ max ]: 

^ r^TY) 1 r) 

+ = tG(0,T] and E G i^max), 

' m{0,E) = mo, E e [E^in, E^ax], 

^m'(f,Emi„) = m'{t,E^^^) = 0, t e [0,T], 

where mo is the initial density. This equation is also called the Fokker-Planck equation in 
physical literature. The detailed derivation of this equation can be found in [14], [29], and 
[30]. 

2.2 The revenues of agents 

Every agent in our model should choose the best emission reduction level r to maximize 
its own net revenues, which consists of three parts, namely, the production revenue, the cost 
of emission reduction, and the net revenue from the emission permits trading scheme. 

Firstly, each player’s revenue arising from the production can be represented by an 
increasing concave function R{Q{t)). Following [20] and [31], we assume that the relationship 
between production and emissions is linear, and the production revenue function can be 
expressed by the following quadratic functional form in terms of emissions: 

, , AE- \E^ 

R(E) = -, 

Cl C2m{t, E) 

where A = Emax, A, and C 2 are constants. Under this assumption, the marginal revenues 
are positive and decreasing. In addition, the production revenue is decreasing with respect 
to the density of agents m, which can be explained by the economical concept “negative 
externality”, that is to say, an agent should face more intense competition and receive fewer 
revenues if it chooses the same state as others’ ones. 

Secondly, the cost of emission reduction should be 

C(.) ^ 

This quadratic form guarantees increasing marginal mitigation cost. 

Finally, the gains from emission permits trading are expressed by 

T{E) = S{t)*{E-Eo), 
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where S(t) denotes the emission permits price and Eq is the initial quota. The trading 
volume of emission permits E — Eq > 0 means that an agent purchases the emission permits 
from the market, and E — Eq < 0 means that an agent sells the unused emission permits to 
others, respectively. 


2.3 The maximization problem and the optimal conditions 


Now, we are in the position to dehne our maximization problem for each agent. That 
is, the objective functional and the constraint condition are as follows: 


maxE< 

n 




AE - \E^ 


r, 


/ —-S(t)(E-Eo) 

Cl + C2m{t, E) 2 

subject to dEt = —Ttdt + adW + dNt{Et), E{0) = Eq, 

where r is the risk-free discount rate, and t = 0 is the initial time. 
According to (2), this problem can be reformulated as follows: 


dt 


(3) 


max 

n 


-^—rt 


0 


AE - lE^ 

Cl + C2m{t, E) 


- S{t){E - Eq) m{t, E)dE 


\dt, 






(4) 


with the initial condition m{0,E) = tuq. 

The solution of the above problem should correspond to the equilibria of mean held 
game, in which every producer is atomized and has rational expectations. 

Next, we will show the process of obtaining the optimal conditions for problem (4). To 
begin with, we multiply v on both sides of equation (2) and integrate by parts to obtain the 
following weak form: 


/ {v{T,E)m{T,E)-v{^,E)m{{),E))dE = 

^ -^min 

for every v G T^^ax] X [0,T]), where 


f^max / Qy 

W7 + 


1 2 dv\ , ^ , 

— T—— 1 mdEdt 


\dt 2 dE^ dE J 


C“(r2) = G C°°(r2); suppw = {x; u{x) 7 ^ 0} C hi [> . 
Then, the Lagrangian of (4) should be 


L(m, r, v) = 


+ 


^—rt 


r*£^max / Qy 


_f 7-2 \ j 

1 


\dt 2 dE'^ dE J 


{v{T, E)m{T, E) — r;(0, E)m{Q, E)) dE. 
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Taking the derivatives of Lagrangian with respect to m, r and n, respectively, we can obtain 
the following system: 


' dv 1 o d'^v dv 


W7 + 


dt 2 dE^ dE 
dv 

dm 1 ^d'^m d 


— Tt— — rv -\r f{E, r, m, t) = 0, n(T, E) = 0, 


(5) 


-c 


{ dt 2 dE^ dE 


+ ^ (-rm) = 0, m(0, E) = mo, m'{t, = 0, 


where f{E,T,m,t) = - S{E - Eq) + ^ • Clearly, we can see that 

this system consists of two coupled equations, one of which is the backward HJB equation 
and the other one is the forward Kolmogorov equation. In fact, the solution of this system 
is the mean held equilibrium of our game. 


3 Numerical methods 

In this section, we will present a numerical method to discretize the hrst equation of 
(5) for the reason that it is difficult to solve the equation analytically. In fact, here a htted 
hnite volume method will be employed. Also, it will be shown that the system matrix of 
the resulting discrete equations is an M-matrix, which guarantees that the discretization 
is monotonic and the discrete maximum principle is satished, such that the scheme has a 
unique solution. Besides, a two-level implicit time-stepping method is used to implement 
the time-discretization. Since the structure of Kolmogorov equation is similar to the HJB 
equation, here we only discuss the latter to save the space. 


3.1 The fitted finite volume method for spatial discretization 


A dehned mesh for is signihcant in the process of discretization. We hrst 

divide the intervals / = (i^min, -^max) into N sub-intervals: 

h '■= {Ei, Ei^i), i = 0,1, • • • , A^ — 1,, 


in which 

For each i 


-^min Eq Ei < • • • < E]\f F/jj^^x- 
0,1, • • • , A^ — 1, we dehne another partition of I by letting 


E,.i = 
* 2 


H-1 


E, 




Ei -I- E, 


i+l 


To keep completeness, we also dehne E_i = i^min and Ej^_^_i = i?max- The step is dehned by 
h = Eij_i — E^_i for each i = 0,1, • • • ,N. 
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Then, for the purpose of formulating hnite volume scheme, we write the hrst equation 
of (5) in the following divergence form: 


dv d f dv ^ , 


( 6 ) 


where 


a = -cr^, b = —r, 

2 ’ 

db dr 

c = r + -— = r — ——. 
dE dE 


It follows from integrating equation (6) over {E^_i, E^_^_i) and applying the mid-point 
quadrature rule to the resulting equation that 


dvj j 

dt * 






Piv)\l 




for i = 1, 2, ■ ■ ■ , — 1, where k = — E-_i, Ci = 

and p{v) is a flux associated with v dehned by 


“t“ fi^i 0 

c{Ei,t), Vi = v{Ei,t), fi 


(7) 

/(Ei,r,m,f), 


p{v) = a-^ + bv. 


( 8 ) 


Now, we focus on deriving an approximation to the flux at mid-point, -Ej+i, of the interval 
li for alH = 0,1, • • • , — 1. Consider the following two-point boundary value problem: 


{av' + = 0, E e /*, (9a) 

v{Ei) = Vi, v{Ei+i) = Vi+i, (9b) 


where 6j_,_i = b{E-_^i,t). A hrst-order ordinary differential equation can be obtained by 
integrating both sides of equation (9a): 


Pi{v) = av' + bi^iv = Cl, 

where Ci is an arbitrary constant and can be determined by solving the above constant 
coefficient two-point boundary problem analytically as follows: 


Pi{v) = Cl = bi^i- 


Ti+i 


e^iEiv- 


oOliE, 


i+1 


■naiEi 


( 10 ) 


where ctj = b^j^i/a. 

By using (10), we can dehne a piecewise constant approximation to p{v) by ph{v) sat¬ 
isfying 

Ph{v) = pi{v) if xeli (11) 

for i = 0,1, • • • , A^ — 1. 
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Thus, substituting (11) into (7), we obtain the following results: 


where 


, L L _ r; 

at 




^ 2 , 2—1 ^i — 


2 ^O-i—iEi _ ^(y.i—\Ei—i 


^Oii—iEi ^cxiEi 

G^ii — b-i -=-=- -|- 6.-_|_i-^^ + Gili^ 

’ ^ 2 1-^i _ — 2 g^i-^i+1 _ ^CtiEi 


Si,i +1 — &i+i 


2 ^ 2 -^^'2 "h 1 ^ 2 ^^ '2 


( 12 ) 

(13) 

(14) 

(15) 


3.2 The implicit difference method for time discretization 

Next we embark on the time-discretization of the system (12). To this purpose, we hrst 
rewrite equation ( 12 ) as 

dv ■ 

—BiV = fik, (16) 


dt 


where 


Di — (ei,i, ei, 2 , 0 , • • • , 0 ), 

Bi (0, * * • ,0, 62 ^ 2 —!, 62^27 0, * * * , 0), t 2, 3, * * * , N 2, 

Dn-1 = (0, • • ■ , 0, e7V-l,Ar-l, CAT-gAf)- 

We select K — 1 points numbered from ti to tx-i between 0 and T, and let T = to, tx = 0 
to form a partition of time T = to > ti > ■ ■ ■ > tx = 0- Then, the full discrete form of 
equation (16) can be obtained by applying the two-level implicit time-stepping method with 
a splitting parameter 9 E [|, 1] to it: 


{eB{E, T, t"+i) + G")n"+' = ef{E, T, t^+y, + (1 - e)f{E, T, t^)l 

+ (G’‘- (I - e)D(E,T,t’‘))v\ 


(17) 


where 


B — (Hi, B2, ■ ■ ■ , Hat-i)^, 

= diag {-h/Atk, ■ ■ ■ , -/jv/Affc)"^, 

for fc = 0,1, ■ ■ ■ ,77 — 1. Note that At^ = t^+i — tk < 0, and denotes the approximation 
of n at t = tk- Particularly, when we set 6 * = |, the scheme (17) becomes the famous 
Crank-Nicolson scheme and is of second-order accuracy; when we set 9 = 1, the scheme (17) 
becomes the backward Euler scheme and is of hrst-order accuracy. 

Given the initial condition of n at t = T, we can solve the values of v at the discrete 
points itk,Ei) by using (17). 

The following theorem declares that the system matrix of system (17) is an M-matrix. 






Theorem 1. For any given A; = 1, 2, • • ■ ,K—1, if | is sufficiently small and c > 0, the 
system matrix of (17) is an M-matrix. 


Proof. First, we note that ejj < 0 for alH, j = 1, 2, ■ ■ ■ ,N — 1, j ^ i, since 

h 




oOiiE-i 


+ 1 


o(y.iEi 


> 0 


(18) 


for any i and any ^ 0. This is because the function e"® is increasing when 6 > 0 and 
decreasing when 6 < 0, where a = ^ and a = Moreover, (18) also holds when —)■ 0. 

Furthermore, from (13)~(15) we know that when c, > 0, for alH = 1, • • • , iV — 1, there holds 


A |(ei,i-i)^^^| + |(ei^j+i) 


fc+i| 


\fc+l I 


+ 


Therefore, D{E,T,t^~^^) is a diagonally dominant with respect to its columns. Hence, from 
the above analysis, we see that for all admissible i, D{E, r, t) is a diagonally dominant matrix 
with positive diagonal elements and non-positive off-diagonal elements. This implies that 
D{E,T,t) is an M-matrix. 

Second, of the system matrix (17) is a diagonal matrix with positive diagonal entries. 
In fact, when \Atk\ is sufficiently small, we have 

H-- > 0, 

—^tk 

which demonstrates that 6D{E, T,t) + G is an M-matrix. □ 

Similarly, we can also discretize the third equation of (5) by using the above method. 
The parameters in discrete scheme can be modified as follows: 

Let b = T, a = b/a and c = 0. Then, the coefficients are given by 


bj 


^Oii — 1 E-i — p 


2 ^Oii — iEi _ ^Cti — iEi — i ’ 


— 1 Ei 


o^iEi 


Cii = b-l— - - -z- - - 

’ ^ 2 ^^i—iEi _ ^(y.i—\Ei—-i 


+ b. 


^”1“ 2 _ ^FX.iEi 


n R. Ghi 


Gi+l ^aiEi+i _ ^diiEi 




These elements can build a matrix D = {Di, D 2 , • • ■ , Dat-i)^, where 


-Di — (ei,i, ei^2, 0, • • • ,0), 

(0, * * * ,0, 62^2—!, 62^27 0, * * * , 0), t 2, 3, * * * , N 2, 

Dn-1 = (0, • • ■ , 0, e7v-i,v-i5 eAr_i,Af)- 


Consequently, the numerical discrete scheme for the Kolmogorov equation reads as 


{eD{E, r, t’^) + = {G’^ - (1 - e)D{E, r, t^+^))m^+C (19) 
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3.3 The algorithm for (5) 


In the above discussions, we have assumed that the control variable r is known, and the 
density m and the adjoint state v can be solved sequentially. However, we can see from the 
second equation of (5) that r is coupled with v. For this reason, we take an iterative method 
to solve the three unknown functions r, m, and v. The algorithm is presented as follows. 


• Algorithm 

Step 1: Give the initial guess of r, and set it as r°. Set a tolerance threshold Tol > 0 
and n = 0] 


Step 2 
Step 3 
Step 4 


Solve equation (19) to obtain m""; 
Solve equation (17) to obtain n”; 
Use v"' to compute Compute 


£” = max||(r*r-(7i‘)”+‘||; 

i^k 


Step 5: If < Tol, let T = m = m"', v = v"', and stop. Otherwise, let n = n + 1, 
and go to Step 2. 


4 Numerical results 

So far, we have been able to show the results of our differential game model numerically. 
We use the following parameter values to solve (5). 

Parameters: T = 1, E m\n = 1, F^max = 5, Ci = 10, C 2 = 0.1, a = 0.3, S = 0.2, Eq = 1, 
r = 0.1. 


4.1 The efficiency of the numerical method 

First of all, we consider the convergence rate of our discretization method to show its 
accuracy and efficiency. Owing to the limitation of space, we only test the adjoint state v. 
Additionally, since the closed-form solution of (5) can not be found, we regard the solution 
obtained by choosing N = K = 512 in both space and time, respectively, as the “exact” 
solution V. We compute the errors in the discrete L°°-norm at the computational hnal time 
step f = 0 on a sequence of meshes with N = K = 2^ for a positive integer n from n = 4 to 
a maximum n = 8. The discrete L°°-norm is dehned as: 

||n^(U,0)-n(U,0)||oo = max \v’"{Ei,0) - v{Ei,0)\, 

l<i<N 

where denotes the numerical solution. The log-log plots of the computed maximum errors, 
along with the linear fitting, are depicted in Figure 1. From the figure we can see that the rate 
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of convergence of in the discrete L°° norm is of the order where h = max {hi). 

Additionally, it demonstrates numerically that our numerical method for (5) governing the 
mean held game is useful and efficient. Some theoretical analysis about the stability and the 
convergence will be discussed in the future works. 



Figure 1: Computed errors in the L°°-norm at t = 0. 


4.2 The solution of the model 

In this subsection, two examples will be addressed to show the equilibrium states under 
different conditions. 

In the hrst one, the initial density of agents mo follows a normal distribution, whose 
mean and variance are 3 and 0.35, respectively. In this case, most of the agents initially 
choose a medium emission level, and the number of agents whose emission level are maximum 
or minimum tends to zero. 

Figure 2 shows the evolution of density for this example. Clearly, we can see from hgure 
2 that as time goes on, the density of agents averagely disperses on the emission interval 
[i^rnin, Fornax] instead of Concentrating at the medium emission level as the initial state did. 
As mentioned above, the influence of mean held composed of all producers on each agent can 
be explained by the economic concept “negative externality”. That is to say, the revenues 
of each agent should decrease if there are a large number of agents at the same emission 
level. To maximize own revenues, each agent tends to choose a different emission level from 
others’. 
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Figure 2: Evolution of m. 

The effect of emission permits price on the equilibrium state is presented in Figure 3. In 
the hgure, the horizontal axis is the emission level E, and the vertical axis is the density of 
agents at the hnal time t = T, namely m{T, E). We plot three results in the cases of S' = 0, 
S = 2, and S = 4, respectively. Obviously, with the increasing of emission permits price, 
the density in lower emission level increases, and the one in higher emission level decreases. 
In other words, more population would like to choose a lower rather than a higher emission 
level when the emission permits are expensive. 

In the second example, we consider a situation in which most of the agents initially 
stay together at a lower emission level. The initial density of agents is piecewise linear on 
[1,2) and [2,5]. In addition, the emission permits price becomes time-dependent instead of 
constant. It is zero at the beginning [0,0.1), then increases from zero to maximum level 
‘S'max = 2 at [0.1, 0.5), and keep this level to the end of the game, which is shown in Figure 
4. 

The equilibrium density of this example is presented in Figure 5. Similar to the hrst 
example, the original obvious aggregation disappears after a period of time. In addition, 
although the initial condition mo is non-differentiable at E = 2, there is no irregular dis¬ 
turbance in the neighbourhood of this point. This can reflect the stability of our algorithm 
more or less. 

Furthermore, Figure 6 shows the effect of the maximum emission permits price S'max on 
the equilibrium. Here three cases, namely S'max = 1, S'max = 2, and S'max = 3, are considered. 
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Figure 3: Effect of S on m(T,E). 


The similar conclusion to the hrst example can be obtained. That is, the number of agents 
in the lower emission level should be increasing as the maximum emission permits price level 
rises. 


5 Concluding remarks 

We present a mean held game model to study the producers’ behaviors in an emission 
permits trading scheme. In our model, there are a continuum of producers, and each producer 
is homogeneous. The equilibrium of this game can be represented by a system containing a 
forward Kolmogorov equation coupled with a backward HJB equation. We then propose a 
so-called htted hnite volume method to discretize the resulted partial differential equations 
and the corresponding system matrix is proved to be an M-matrix. The efficiency and the 
usefulness of this method are illustrated by the numerical experiments. Our results show 
that the convergence rate of our method is nearly 0{h‘^). Besides, due to the externality, 
each agent tends to choose a different emission level from the others’. Finally, we find that 
the number of agents in lower emission level should be increase when the emission permits 
are expensive. 

We anticipate that our methodology from the perspective of partial differential equations 
combined with numerical methods can make a few contributions to the solving of complex 
problems arising from the interdisciplinary helds. 
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Figure 4: Evolution of emission permits price S{t). 



Figure 5: Evolution of m. 
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Figure 6: Effect of S'max on m(T, E). 
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